Methods for the Study of Transverse Momentum Differential Correlations 



Monika Sliarma and Claude A. Pruneau 
Physics and Astronomy Department, Wayne State University, Detroit, MI 48201 USA 

We introduce and compare three differential correlation functions for the study of transverse 
momentum correlation in p + p and A + ^4 collisions. These consist of inclusive, event-wise and a 
differential version of the correlation measure C introduced by Gavin [l|] for experimental study of 
the viscosity per unit entropy of the matter produced in ^4 + ^ collisions. We study the quantitative 
difference between the three observables on the basis of PYTHIA simulations of p + p collisions and 
A + A collisions consisting of an arbitrary superposition of p+p collision events at ^/s —200 GeV. We 
observe that inclusive and event-wise correlation functions are remarkably identical to each other 
where as the observable C differs from the two. We study the robustness and efficiency dependencies 
. . . of these observables based on truncated Taylor expansions in efficiency in p + p collisions and on 

OO ' the basis of Monte Carlo simulation using an adhoc detector efficiency parameterization. We find 

, that all the three observables are essentially independent of detector efficiency. We additionally 

^— ' ■ study the scaling of the correlation measures and find all the observables exhibit an approximate 

l/N dependence of the number of participants (N) in A + A collisions. Finally, we study the 
impact of fiow-like anisotropy on the inclusive correlation function and find fiow imparts azimuthal 
modulations similar to those observed with two-particle densities. 
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I. INTRODUCTION 



Measurements of two- and multi-particle correlations have proven a powerful tool to study heavy ion collision 
dynamics. In particular, measurements of two-particle correlations in the form of two-particle densities as a function 
of the particle pseudorapidities and azimuthal angles in A -\- A collisions have revealed new correlation features not 
found in p + p interactions. Specifically, azimuthal two-particle correlations in central A -I- A collisions exhibit novel 
"away-side" structures, commonly referred as the away-side dip, and a near-side structure called the "ridge" . First 
measurements of the "dip" and "ridge" were reported on the basis of inclusive charged particle correlations [2| . More 

— , recently, measurements were extended to include identified charged and neutral particles, studies of correlations in 
' various momentum ranges and system size dependencies Q . Measurements of event-by-event fluctuation observables 

such as net charge and transverse momentum provide complementary information. Studies of net charge and transverse 
■ momentum fluctuations show that while the magnitude of the fluctuations observed in A + A collisions differs from 
QQ ' that observed in p-\-p collisions (after proper scaling to account for the number of collision participants), no dramatic 
. suppression of the net charge fluctuations is observed that might signal the formation of a Quark Gluon Plasma 
^ ' (QGP) [1, Q. Similarly, only small differences are observed in the magnitude of pt dynamical fluctuations in A A 
• ^ . collisions relative to those observed in p -\- p interactions after scaling by the number of participants and the particle 
' average pt ■ The measured dynamical pt fluctuations in A -I- A thus do not provide evidence for the strongly enhanced 
^ ,, fluctuations expected for nuclear matter near the tri-critical point In contrast, however, the width of the charge 

- - ' balance function exhibits a signiflcant reduction from peripheral to central Au -\- Au collisions that is consistent with 

the predictions of delayed hadronization when a QGP is formed Q . The reduction of the balance function measured 
as a function of pseudorapidity may, however, result from radial flow, resonance decays, or other effects. It is thus 
important to seek additional observables to obtain a better understanding of the A-\- A collision dynamics. 

In this work, we consider measurements of differential transverse momentum two-particle correlations as a function 
of the relative pseudorapidity and azimuthal angles of measured particles. A certain latitude exists in the calculation 
of such correlation functions. Indeed, as for integral correlations @, Qj one can use both inclusive and event-wise def- 
initions. One can also consider dynamical fluctuations as the difference between measured fluctuations and statistical 
fluctuations i.e. those expected for a purely Poisson system. Extensions of fluctuation variables Q and [loj 
are also possible. In this work, we study three differential observables namely, inclusive and event-wise correlation 
functions and a differential version of the measure C, used by Gavin [O, [l^l to estimate the viscosity per unit entropy 
of the matter formed in A-\- A collisions. The observable definitions and notation used in this work arc presented in 
Sect, ini In Sect. IIIIl we focus on the definition of these observables for the study of transverse momentum corre- 
lations between two particles as functions of the relative pseudorapidity, and azimuthal angles of the particles. The 
goal of these observables is essentially to measure the covariance of two particle momenta by averaging the product 
of the differences between these momenta and the mean momentum. Given their different definition, one expects, 



2 



in practice, they should measure different numerical values. The quantitative differences between these observables 
are assessed in Sect. IIVI on the basis of generic arguments and PYTHIA simulations. Section |V] presents predictions 
based on the event generator PYTHIA Sect. I VII discusses the dependence of these observables on the number of 
participants in ^ + A collisions. Elliptic flow effects on the magnitude and azimuthal dependence of these observables 
is discussed in Sect. I VIII The experimental robustness of these observables is studied in Sect. IVIIII This work is 
summarized in Sect. IIXI 



II. NOTATION AND DEFINITIONS 



We consider measurements of one- and two-particle production cross-sections. We are specifically interested in the 
measurements of two-particle correlations. Let and yi represent arbitrary observables for two particles, z=l, 2. 
These observables may be the momenta, azimuthal angle, pseudorapidity etc., of the particles. Let pi (xi) and 
P2 {xi,X2) be one- and two-particle densities, respectively, defined as functions of these observables. 

Pi ix,y) = -j—r- (1) 
axdy 

dN2 

P2ixi,X2,yi,y2) = , , , , (2) 
axiax2dyiay2 

Discussions in this paper are limited to two-particle correlations but can be straightforwardly extended to higher 
order multi-particle correlations. For simplicity, and illustrative purposes, assume the particles are indistinguishable. 
Integrals of the above densities over the space spanned by observables (xi,yi) yield the number of particles and pairs 
of particles, respectively, in the domain of integration. 

pi {x, y) dxdy = / —^dxdy = N (3) 
J dxdy 

f f dN2 

/ P2{xi,X2,yi,y2)dxidx2dyidy2^ / 3 — -, — -. — —dxidx2dyidy2 = N {N - \) (4) 
J J dxidx2dyidy2 

where N is the average number of particles produced in the region of acceptance of interest. 

Experimentally, the number of particles in a given acceptance fluctuates event-by-event owing to the stochastic 
nature of the production processes. The one- and two-particle densities may then be viewed as an average of the 
number of particles taken over all measured events. We label events with an index a = 1, ...,N^y, where Nf,.^ is the 
number of events in the data sample. One measures the number of particles, na{y) at a given location y in a bin 
of width Ay for each event a. The one- particle and two- particle densities are then obtained by averaging over the 
event ensemble, i.e. the total number of events N^v 

N 

/\ 2. 1 " " 

Piix,y)dx —{n{y)) = ——Y^^Ua^iy) (5) 

N 

P2 {xi,X2,yi,y2)dxidx2 = {n{yi)n{y2)) = -^Jf- X] '^^^^i^'^^^^^) (6) 



where (O) represents the event-ensemble average of O. 

In this work we describe measurements of fluctuations and a differential correlation observable. This requires 
measurements of covariances (Aa;iAa;2) where Ax = x — {x) and (x) is the mean value of x, which is a function of 
coordinates yi and 1/2 ■ Such a function of yi and 1/2 corresponds to a Aa:iAa;2 average taken with respect to the 
two-particle density given by Eq. 2. To denote such averages we use the notation pn''^'' (y), where g(x) is some function 
of the coordinates xi, X2 and n = 1,2 signifies the one- or two- particle density, respectively. For example, the average 
value of a; as a function of y, as defined below, denoted pfiy), is defined and calculated as: 

= !p.ix^y)dx - = ^^r-— 
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where Xa.i{y) is the value of x for particle i at position y in an event a. similarly the two-particle covariance (Aa:iAx2), 
a function of j/i and ?/2 is denoted P^^'^^^^{yi,y2) and is defined and calculated as: 



A2;iAx2 



iyuy2) 



{Axiyi)Axiy2)) 



^ J P2 {xi,X2,yi,y2) A.xiAx2dxidx2 
J P2ixi,X2,yi,y2)dxidx2 

E E E (xaAyi) - {x{yi))){xa,j{y2) - {xiy2))) 
na{yi)na{y2) 



(8) 



The number of pairs in an event a is given by the product of the number of particles, na{yi)na{y2), in bins of 
width Ay at positions yi and 7/2- However, the number of pairs is reduced to na{yi) {na{yi) — 1) for identical bins 
yi = y2 (indistinguishable particles). 

The above expressions correspond to what is commonly referred to as inclusive averages. An alternative averaging 
method is also available which corresponds to event-wise average, event-wise averages of x and Aj;iAa;2 are herein 
noted as pf{y) and P2^^^^^ iyi,y2), respectively. They are calculated as follows: 



p: (y) 



— Y — 



Nev ^ «a(yi) na{y2 



(9) 
(10) 



where 



i=l 



{x{yi))) {xa,j{y2) - {xiy2))) 



(11) 



For integral correlations where the above averages are taken over a single wide y bin, the above expression (Eq.lO) 
reduces to fluctuation variables reported in recent works 0, [3] . 

Gavin suggested recently [T^ that the width of transverse momentum fluctuation distributions measured as a 
function of pseudorapidity may be used to estimate the viscosity of the matter formed in A + A collisions studied at 
RHIC and LHC. 

Here we consider a differential version of the correlation C, and write 



C = 



E E E Xa,t{yi)Xa,j iy2) 
a—1 i—1 i^j — 1 



E E Xa,r (2/1) ' 
a—1 i—1 



X; ria {yi)na (^2) 



V 



E (2/1) 



/Af„"»(y2) \ 

E E Xa.j (2/2) > 
a=l j=l 



J 



V 



E n-a iy2) 



(12) 



A number of other observables have been used recently to carry out measurements of fluctuations and are well 
documented in the hterature. The quantity $pj, defined in was used by CERN experiments NA49 and CERES 
to measure transverse momentum fluctuations in Pb + Pb collisions. The quantity a^^ was used by STAR to 
measure fluctuations in Au + Au collisions ■ There is an approximate equivalence between these different integral 
observables as shown in UM 



III. TRANSVERSE MOMENTUM DIFFERENTIAL CORRELATIONS 



In this section, we focus our discussion on integral and differential transverse momentum correlations. Integral 
correlations, (Apt Apt), may in principle provide an estimate of event-by-cvcnt temperature fluctuations which, in 
turn, may provide means to determine the heat capacity of the medium formed in A + A collisions (itI . [isl . [l^. 
Additionally, one can also use the measure C, defined in Sect. 2, to study the transverse momentum correlations 
0, [m, • Wc argue that a differential version of all the observables, as defined in the previous section, may provide 
an equivalent or even better way to estimate the medium viscosity. However, calculation of viscosity using these 
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observables is not presented in this work. Clearly, differential correlation functions provide more information than 
integrals over a wide range of pseudorapidity. Therefore they may enable a deeper understanding of the collision 
dynamics at play in A + A collisions. As a specific example one notes the discovery of the "ridge" in two-particle 
(density) correlations measured in A + A collisions. It is thus interesting to consider differential transverse momentum 
correlations as a function of differences in particle pseudorapidity and azimuthal angles (denoted as A77 and Aip, 
respectively). 

To simplify the notation, we use p rather than the "usual" pt to denote transverse momentum. Wc define the 
transverse momentum correlations of interest by substituting p for x in Eq. [7] through Eq. 1121 and include a 
dependency on the particles' pseudorapidity, 77, and azimuthal angle, ip. The inclusive momentum correlation function 
is defined as 



gApiAp2 = 1 P-^ {Pi,P2, At?, A(p) ApiAp2dpidp2 



J p2 {pi,P2, A?/, A(p) dpidp2 



(13) 



and calculated using the following expression: 



ApiAp2 



(A?/, A</7) 



JVci. nnC'U.Vl) "0(172, V2) 

EE E (Pa,*(»7l,</'l) - (p(m,Vl))) (Paj('72,</'2) - (P(?72, V2))) 

a— 1 i—1 jz^i—l 

a=l 



(14) 



where one applies the conditions A?/ = ?yi — 7?2 and Aip — ipi — Lp2. na{rii,^pi) is the number of particles detected 
in a given event a at pseudorapidity rji and laboratory angle ipi. Pa.iiViT ^i) represents the transverse momentum of 
particle i from event a. {p{rji, ipi)) is the inclusive average particle transverse momentum at rji and fi. 



E E PaAV.'P) 

{p iv, </')) = Pi iv, v>) = jv'r^ 

E naijq.ip) 



(15) 



We include dependencies of the average transverse momentum, on the laboratory angles Lpi and ^p2 to account for 
possible instrumental effects. The mean p is independent of azimuthal angles Lpi and ip2 for unpolarized beams. The 
detector acceptance and particle detection efSciency may, however, vary with these angles and must thus be accounted 
for in practice. 

The event-wise transverse momentum correlation is calculated as follows. 



-p^^^"^^^ {Ar^,Ap) = ^Y. 



1 



'ev ria imiVl) ria (?72, V2) 

rac(r;i,Vl) nc.('/2,¥'2) 

X 51 X! iPo',iiVl,Vl)- P^liV,f))iPa,j{V2,(P2)- P^iiV,^)) 

i=l i#i=l 



(16) 



and uses the event-wise average 



(17) 



Similarly C can be written as: 



C 



EE E Pa,i {111,^1) Pa,] {m,V2) 



a—1 i—1 



E ria ivi,Vi) n-a (f?2,'P2) 



/ Wet, riair]i,ipi) 

E E Pa,z im,y^i) 

a=l i=l 

E ivi,'Pi) 

a=l 



E E Pa,3 {m,V2) » 

E "^a (V2,'P2) 
a=l 



(18) 
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IV. COMPARISON OF DIFFERENTIAL CORRELATION OBSERVABLES 



The inclusive and event-wise differential correlation observables, as well as the correlation measure C, are designed 
to study the magnitude of correlation in terms of some specific kinematic variable x, where x may be, e.g., the particle 
transverse momentum or pseudorapidity. The three observables represent averages of the difference between the x value 
of each particle relative to the mean value of x. The three observables, however, employ distinct averaging methods 
and in general lead to different numerical results. None of these averaging methods can be considered "superior" 
in any sense. However, due to practical or technical reasons, studies of correlations are typically reported u sing one 
type of averaging only. Still other averaging methods were advocated and used in recent works jl, [lo, [Ml, [13, [U [11] . 
Experimenters, therefore, have a choice of observables and have in fact reported data using many of these averaging 
methods. The different methods used by several experiments produce results that cannot be trivially compared from 
one experiment to another. The issue is compounded by the fact that these observables have different dependencies 
on the detection efficiency and acceptance. It is therefore of considerable interest to determine how different are the 
numerical values produced by these different averaging methods and correlation observable definitions. One can also 
find out whether values measured by an experiment using one observable can be meaningfully compared to those of 
different experiments or measurements using other observables. 

We first note that if particle production follows perfect Poisson statistics, the correlation functions defined in the 
previous section are, by construction, null. The inclusive, event-wise and C observables are thus trivially identical 
in the Poisson fluctuation limit. However, due to energy, momentum and quantum number conservation, particle 
production in general is a non-poissonian stochastic phenomenon, i.e. the correlation functions are non-zero. Below 
we investigate how different they can be in practice. 

In order to investigate the differences between the inclusive, event-wise and C observables, we consider a decompo- 
sition of the one- and two-particle densities using fixed multiplicity densities. We write 

oo 

PliP,V) = Pi^)Phm{P,V) (19) 

m=0 

oo 

P2{pi,m,P2,V2) ^ X! P{fni,m2)p2 , rni, 7712 (pi, 1)1 iP2, 1)2) 

(20) 

7711 ,7712—0 

where P[m) expresses the probability of having a multiplicity 'm' in the measured 77 and Lp bin. Similarly, P(mi, 7712) 
corresponds to the probability of finding mi and 7712 particles simultaneously in bins (771, ipi) and (772, ^2)1 respectively. 
Pl,m{Pi'n) and P2,m {pi,'ili,P2,ri2) arc the single and pair densities for a fixed value of multiplicity m (7711,7712). Note 
that in the above expressions we have included a dependency on ?y only for brevity ((p dependencies can obviously 
also be included). By construction, one has 



Pi.niip, r])dpdr] = m (21) 



P2.mi,m2{Pi,miP2,'r]2)dpidp2dr]idrj2 = mim2 (22) 
Integration restricted to particle momenta yields the pseudorapidity densities: 

Pi,m{PiV)dp = Pi,r,i{v) (23) 

{Pl,Vl,P2,V2)dpidp2 P2, 7711 ,7772 {Vi,V2) (24) 

The averages of the transverse momentum, p, and ApiAp2 at fixed multiplicity 'tti' are 

/ (,) ^ IpiMp.v)pdp ^25) 

Pl..m[V) 

pApAp r ^ ^ / ^2.7771,7772 iPl,m,P2, 772 ) Agl Ap2 ^Pl ^^2 ,26) 

^' ^ P2,777i,m2(m,»?2) 
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Given the multiplicity m (mi, 7712) is (are) fixed, the above expressions are valid and identical for the inclusive and 
event-wise averaging methods. Thus for fixed multiplicity densities the inclusive and event-wise averages are truely 
comparable for analyzing actual densities such as occur in p-\-p and ^ + A collisions. Based on Eqs. [T3 l fT4 l [T9l fc [20 l 
the inclusive pt and Apt Apt means are expressed as follows in terms of the fixed multiplicity densities: 

P'liv) = (27) 

m=0 
oc 

E P(mi,TO2) / P2,mi,rn2iPl,Vl,P2,ri2)^Pl^P2dpidp2 
ApiAp2/ N _ mi,m2=0 



J2 P{fni,m2)p2,ra{r}iim) 

mi,m2— 

00 

E -P(mi , TO2 )p2,mi ,m2 [m i ^2 )pf'^'f'^^ [Vl , m) 
mi,m2— 

oo 

rni ,m2— 



(28) 



The mean pt and the /^pt^Pt fluctuations within a single wide range of pseudorapidity are given by 

oo 

E ^(™)"^pf „(f?) 

(P> = ^ (29) 

E P{m)pi,^{vi) 

m=0 
oo 

E F(m)m(m- l)pfP^P(7]i, 7^2) 

(ApAp) = ^ ^ (30) 

^ P{m)m 

m— 

The event-wise mean pt and the AptApt arc calculated in a similar fashion. 

??i \ oc 

i=l I m=0 

P2''^''i.m,m) = ( E E Ap,(77i)Ap,(,?2)) 



I 2=1 ]^t=l 



E ■P("^i''^2)p^^'^^„,(?7i,'?2) (32) 

mi ,m2— 



For a wide bin of integrated pseudorapidity this becomes 

I -j^ ra \ oo 

-EpO = E^("^)(p)™ (33) 
^—[^ E P^^'^- ) = E ^("^) (^P^?')™ (34) 



m 
i i=i 



ApAp 



m(m — 1) 



Lastly, C expressed in terms of fixed multiplicity densities is calculated as follows: 

^ ^ {mjm^l) {piP2)J _ ^^^2 



E P(m)?7l(77l- 1) {piP2)m 

(p)' (35) 
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Though the inclusive, event-wise and C quantities are manifestly different, they exhibit interesting similarities. We 
first consider the difference between the inclusive mean, (p) , and the event-wise mean, pt , for large integrated ranges 
of pseudorapidity. Assuming (p)^ is roughly constant over the multiplicity range of interest, one finds 



m 



oc oc 
m— m— 

i.e. the inclusive and event-wise mean pt should be approximately equal. This approximation is most likely valid 
for large domains of integration in pseudorapidity and in central A + A collisions, which yield large multiplicities. In 
this case, it is reasonable to expect that (p),„ varies slowly with m and one thus expects the difference between the 
two observables to be the smallest. On the other hand for small domains of integration or small multiplicity (p).^ 
is more likely to depend on m and consequently one expects the two average methods to exhibit larger numerical 
differences. 

The inclusive, event-wise and C two-particle correlations can be approximated in the same manner. One finds for 
integral, or wide bin averages: 

oc 

^ P{m)m{m — I) 

{Ap,Ap,) « (ApiAp^)™ = (ApiAp2)™ (38) 

OO 

ApiAp2 « (ApiAp2),„ E ^("^^ {^Pl^P2)m (39) 



E P{m)m{m - 1) (ApiAps)™ + (p) 



\ 2 



= (P) 

Im) 



2 



(ApiAp,) - ) " +^ — , ip) 

m m 



2 



1 (m(m — 1)) , ,9 / (7Ti(m — I)) \ 
« (ApiAp2) ^ ) " + {pf ^ \ . " - 1 (40) 
[m) \ {m) j 

where we have again assumed both and (ApiAp2),„ vary slowly with 'm'. One finds that in the limit where 
(p)^ and (ApiAp2)m are independent of m, the inclusive and event-wise observables are strictly equal. In practice, 
the equality is only approximate given the average particle momentum is a function of event multiplicity. Additionally, 
given that for large values of m, and relatively narrow range of multiplicity values, the ratio is, typically near 

unity. Therefore, one concludes that C is also roughly equal to (Api Ap2) for large multiplicities. The above arguments 
and approximate equalities hold for integral correlations (or fluctuations) but are weaker for differential correlations. 
Indeed, in the case of differential correlations, the number of particle in a given pseudorapidity bin is quite small. The 
approximations used in Eqs. [2T]&[H]are thus likely to be inappropriate and one should expect significant differences 
between the three observables; particularly between the inclusive average and C. Any comparison of results from 
experiments reported on the basis of different observables is thus at best qualitative. We explore these differences 
further in Sect. IVlon the basis of PYTHIA simulations. 



V. MODEL PREDICTIONS 

We use PYTHIA (version 6.22) [T^ to simulate p+p colUsions at = 200 GeV and carry out calculations of trans- 
verse momentum correlations based on the observables introduced in Sect IIIII Five million minimum bias PYTHIA 
events were integrated to produce plots shown in the following. Figures [IJa-c) show the inclusive, p^P^^P^ (^Ari,Aip), 

event-wise, p^^'-'^^^ {Arj, Aip) and C correlation functions, respectively, for particles in the pseudorapidity range 
\ri\ < 1 and transverse momentum range 0.2 < pt < 2.0 GeV/c. The two-particle correlations are plotted as a function 
of the particles' relative pseudorapidity, Ar], and azimuthal angles, Aip, using 31 and 36 bins, respectively. 
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The three correlation functions show quaUtative similarities: all three observables exhibit a peak structure near 
Ary = A(p = (near side) and a ridge-like structure at A(p = tt (away side) which extends over the full At; range. 
These shapes are qualitatively similar to that of density correlations predicted by PYTHIA While the inclusive 
and event-wise correlation functions are identical to each other with a difference of only 0.2% (statistical errors are of 
the order of ±0.001), the observable C exhibits characteristically different strength and shape. In contrast to inclusive 
and event-wise observables, the away side ridge-like structure has larger magnitude than the near side. One also finds 
that C has a strength five times larger than the other two observables. A detailed study of C reveals that the shape is 
largely determined by density correlations, ^^^^^ . In contrast to inclusive and event-wise measures, which are designed 

to minimize the number density contribution C is sensitive to the variation of the pt of the particles as well as the 
number density. As a result we conclude that the differences between the three observables stem from their specific 
definitions and hence they lead to very different numerical results. 




FIG. 1: (Color online) Comparison of (a) inclusive, ''^ (A?;, A<^), (b) event-wise, ''^ ''^ {Ar],A(p) and (c) C transverse 
momentum correlation functions obtained with p -\- p collisions at -y/s — 200 GeV generated with PYTHIA. 

Measurements of two-particle correlation functions, p2(A77, A(^), in Au -\- Au collisions reveal the presence of a 
strong ridge-like structure on the near-side. A(p w 0, and small correlation yield on the away-side, Aip w tt, while 
ellip tic flow appears as a rather modest contribution to the overall yield of the correlation function. Authors of Ref. 
[2a | have suggested the change in the shape of the correlation functions, in which a ridge-like structure shifts from 
Aip « TT to Aif « 0, might in part be caused by strong radial flow present in mid to central A -\- A collisions. It 
is thus interesting to test the radial flow scenario by considering the effects of radial flow on transverse momentum 
correlations. Figures [2fa-c) show the inclusive, event-wise and C correlation functions calculated for radially boosted 
p -\- p events. All particles produced by a given p -\- p event are boosted radially in the transverse plane with factor 
/3 = 0.3. Radial flow imparts kinematic focussing, as a result of which all the particles from each p + p collisions 
are pushed in the same direction and hence become correlated in azimuth. Indeed one finds a significant correlated 
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yield at large pair separation in pseudorapidity and narrow separation in azimuth on the near side of all the three 
observables. The inclusive and event-wise measures are qualitatively similar with a peak and ridge-like structure 
on the near side and diminished correlations on the away side. However, C differs quantitatively. The away side 
correlations are reduced significantly and the near side peak-like structure is more prominent vis-a-vis inclusive and 
event-wise observables. Also, the strength of C for boosted p-\-p events is nearly 10 times larger than the strength of 
the other two observables. These simulations show that radial flow in A + A collisions can produce a near-side ridge 
extending in pseudorapidity in transverse momentum correlations as well as in two-particle density correlations. 





FIG. 2: (Color online) Comparison of (a) inclusive, p^^^ (A77, A<^), (b) event-wise, p^^^ ^'^ (Ary, A(^) and (c) C transverse 
momentum correlation functions obtained with radially boosted (/3 = 0.3) p + p collisions at = 200 GeV. 



VI. SCALING WITH NUMBER OF PARTICIPANTS IN HEAVY ION COLLISIONS 

In this section, we consider the Api Ap2 dependence on the number of participants in heavy ion collisions. We use as 
a reference collisions where nucleon-nucleon interactions are independent of one another and we neglect re-scattering 
of the produced particles. 

Under such assumptions, the two-particle density in A + A collisions at fixed number of participants. N, may be 
written as follows: 

P2.AA = N{N - l)pi,ppPl,pp -I- Np2,pp (41) 

The inclusive differential transverse momentum correlation function at fixed N in A+A collisions. P2^2a^^ ivi: V2) 
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is then: 



'2,AA ^'^I'^sj^ N{N~l)pi,pp{r^i)pi,pp{'n2)+Np2,pp{m,V2) 
^ Np2^ppiVi,V2)pt:fP^{Vi,V2) 

N{N- i)pi.pp{m)pi,ppim) + Np2.^pp{m,V2) 

where we have included, for simplicity, only dependencies on the pseudorapidity of the particles. The simplification 
obtained on the second line in Eq. [42] arises because the average of Apt (first term of the numerator) is by construction 
null. Experimentally, one cannot constrain the collision impact parameter, so one must account for fluctuations in 
the number of participants. The inclusive differential correlation function is then 

Ap,A,., , {N)p2^pp{m,m)p^,^f''^im,m) 



' ' {N{N~l))p,,ppir„)p,,pp{r,2) + {N)p2,pp{m,V2) 



{N{N-l))p^,pp{r1,)p,,pp{ln)|p2,pp[m.m) + {Nf^■^^ ^'^^''^'^ ^ ^ 

Fluctuations in the number of p + p interactions have Poisson statistics in the context of our independent collisions 
model. It implies: {N{N - 1)) = {Nf > {N) Additionally, since p^^P\m,m) / P^^Km)p^^\m) ~ 1 one expects 

P^'^^''- im .m)^ -^fZt'' ' (44) 

i. e. the inclusive transverse momentum correlations m A -\- A collisions should be approximately proportional to 
those measured in p + p collisions and inversely proportional to the number of participants. 

We check this result by artificially constructing + ^ events on the basis of a superposition of independent 15 
p+p collisions generated with PYTHIA. Figures[3l^a-c) display correlation functions obtained with A^= 15 independent 
p + p collisions modeling A + A collisions with 30 participants. A total of two million minimum bias events were 
integrated to produce these results. 

The most prominent feature observed for all the observables is the dilution of signal strength. The signal strength 
is diluted by a factor of 15 as compared to that found in p + p collisions (shown in Figure [1]). This effect is attributed 
to 15 times larger multiplicity produced for 15 independent p+p collisions. The shape of inclusive and event-wise 
observables, with a peak-like structure on the near side and a ridge-like structure on the away side is still observed, 
as previously observed for Figure [IJa,b). Similarly we observe the same shape of C, i.e. larger amplitude of the away 
side ridge vis-a-vis the near side peak. 

We also present results on ^ -f- A collisions produced by the superposition of = 100 p + p collisions, i.e., 200 
participants. Two million minimum bias events were integrated to present the results in Figure [Ja-c). As expected, 
we observe the signal strength is diluted by a factor of 100. 

We next consider the effect of radial flow on artificially generated collisions. As for Figures [5ja-c), we produce 

A + A collisions consisting of N = 15 independent p + p interactions. Particles of a given interaction are now boosted 
radially in the transverse plane with (3 = 0.3. The dominant features remains the same as discussed previously for 
Figures [2{a-c), i.e., observation of diminished away-side ridge-like correlations and enhanced near-side correlation for 
all the three observables. Therefore, on the basis of simulated p + p and A + A collisions, we conclude that strong 
radial flow should shift the away side structure to the near side in A + A collisions and thereby produce a ridge in 
momentum correlations. 



VII. EFFECT OF FLOW ON TRANSVERSE MOMENTUM CORRELATIONS IN HEAVY ION 

COLLISIONS 

Flow plays an important role in heavy ion collisions. It particularly affects particle correlations in azimuth. The 
effect of elliptic flow on azimuthal (density) correlations is well documented [2^. Here, we consider its effect on trans- 
verse momentum correlations studied as a function of relative azimuthal angles and pseudorapidity. The calculation 
method is similar to that used in 

We model the probability density function (pdf) of particle emission as a function of pseudorapidity, azimuthal 
direction and transverse momentum, relative to reaction plane as follows. 

Fi (r?, p| *) = Pi (r/, p) 1 1 + 2 ^ w„ (r/, p) cos (n {if - *))| (45) 
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FIG. 3: (Color online) Comparison of (a) inclusive, p^fi^P^ (A77, A</p), (b) event-wise, p^P^^P^ {Ar],A(p) and (c) C transverse 
momentum correlation functions obtained with superposition of A'' = 15 independent p+p collisions at ^/s — 200 GeV generated 
with PYTHIA. 



^ represents the reaction plane azimuthal direction. 77, tp, and pt are the pseudorapidity, azimuthal emission angle 
and transverse momentum of the particle, respectively. Pi{ri,p) is the pdf of the particle at a given pseudorapidity 
and transverse momentum. Explicit knowledge of this pdf is not required in the following. Assuming the reaction 
plane is not measured, the two-particle pdf is given by 

-P2(?/l,<^l,Pl,'72,</'2,P2) = Pl{in^Pl)Pl{m^P2) X |l + 2^W„(77i,pi)w„(772,P2) COS (n {ipi - ^Pt))^ (46) 

We use Eq. [T3]to calculate the inclusive ApiAp2 correlations: 

2EK('?i) - {pim))vnim)) - {p{m))vnim))cos{n{ipi - (p2)) 
p^p^'^p^(T]i,Lpi,r]2,Lp2) = — . , „^ — r~\ — r~\ — n v; (4"^) 

1 + 2}2Vnim)Vn{ll2) cos [n [ipi - 1^2)) 
n 

where Vn {rj) and wJJ' (77) are average and pt weighted average of flow coefficients function as a function of pseudora- 
pidity, respectively. 

Vniv) = / Pn{V.P)Vn{-n,p)dp (48) 

1 



Univ) = „ \ \ I Pn{v,P)vn{r],P)pdp (49) 
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FIG. 4: (Color online) Comparison of (a) inclusive, p^P^^P^ {AT],A(p), (b) event-wise, p^Pi^^a (Ary, A<^) and (c) C transverse 
momentum correlation functions obtained with radially boosted (/3 = 0.3) A'^= 100 p + p collisions at ^/s = 200 GeV, generated 
with PYTHIA. 



and 

We find the magnitude of the Apt Apt correlations associated with a flowing medium are determined by the flow 
coefficients weighted by the transverse momentum and the averaged transverse momentum of the system. The 
correlations exhibit cosine modulations in azimuth. It is observed in Au + Au collisions at RHIC that flow coefficients 
and average transverse momentum arc both functions of the pscudorapidity of the particles [28l . [29| . We therefore 
expect, based on Eq. [TH the Apt Apt two-particle correlations associated with flow are also a function of the particles 
pseudorapidity difference. Azimuthal anisotropy coefficients measured at RHIC are largely dominated by elliptic flow 
(W2) coefficients. Flow modulations of the inclusive correlations should then also dominated by second order harmonic 
coefficients. 



VIII. EXPERIMENTAL ROBUSTNESS OF CORRELATION FUNCTIONS 



Measurements of correlation observables in p -f- p and A + A collisions require one applies corrections for the 
limited acceptance and detection efficiency. In some cases, correlation functions may be formulated in such a way 
that detection efficiency partly cancels. This was discussed already by a number of authors for various correlation 
observables. See for instance Ref. [lOl for a discussion of the robustness of factorial moments and the observable, 
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FIG. 5: (Color online) Comparison of (a) inclusive, p^^^^^^ {AT],Aip), (b) event-wise, p^P^^P^ (Ary, A(p) and (c) C transverse 
momentum correlation functions obtained with radially boosted (/3 = 0.3) A'^= 15 p + p collisions at ^/s — 200 GeV, generated 
with PYTHIA. 



_dym used in measurements of charge correlations. In this section, we investigate the extent to which the transverse 

momentum differential correlation functions, defined in Sect. IIIIl are robust observables. We use simple simulations 
to explicitly verify the robustness of these observables in contexts where efficiencies are non-trivial functions of the 
particle kinematics. 

Neglecting instrumental effects such as track splitting, ghost tracks etc., the single-particle and pair yields, Ni(x) 
and N2{xi,Xj) measured at some pt, rj, Lp, are the products of the one- and two-particle densities by efficiencies ei{x) 
and e2{xi,Xj) for measuring one- and two- particles, respectively. 



Ni{x) = ei{x)pi{x) (51) 

N2{,Xi,Xj) = £2(Xi,Xj)p2{Xi,Xj) (52) 

In the above equations, x represents the measured variables pt, r\, ip. The 'raw' momentum correlation, 
/^2^measured('?i' '^s), is thus a fuuctiou of the detection efficiency. We show in the Appendix that particle losses 
due to limited efficiency in general lead to a modification of the fixed multiplicity momentum averages. Measure- 
ments of average momentum are consequently intrinsically non-robust. However, one also finds that if the average 
Pt is independent of, or varies slowly with, the event multiplicity, one can apply Eq. [52] to evaluate the correlation 
functions. 

ApiAp2/ N J £2{m,'Pl,Pl,V2,'P2,P2)P2{Vl,'Pl,Pl,V2,'P2,P2)PlP2dpidp2 

P2,raw {Vl,Vl,m,V2) = f 7 r^—, (53) 

J £2{m,^l,Pl,m,V2,P2)P2{Vl,'fl,Pl,m,'i^2,P2)apidp2 
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Obviously, if the efRciency is independent of the measurement coordinates, i.e. e2ir]i,Pi, V^ij ^2,P2, </'2) = constant, 
then efficiencies cancel out in Eq. 1531 the inclusive transverse momentum correlation is thus a "robust" quantity. In 
general, however, the efficiency is a function of the measurement coordinates. This implies corrections are required 
to account for detection efficiency. We note that Eq. [52] simplifies considerably if the efficiency can be factorized 
into separate functions of momentum and pseudorapidity. It reduces to p^^^'^^^ {i11jV2), *-e- is robust to first order, 
if the efficiency is uniform across the pt acceptance of the measurement. Large detectors at RHIC and LHC have 
detection efficiencies that are reasonably uniform across their acceptance. However, various instrumental effects, such 
as detector boundaries or defective components, introduce small non-uniformities. We model these as follows: 



£1(2:1) = ei,o(l + ^i(a;i)) (54) 
£2ixi,X2) = £2,0 (1 + 52(2:1, 2:2)) (55) 

where ei^o a-nd £2,0 correspond to the average efhciencies for detecting one- and two-particles, respectively. The 
functions di{xi) and 62{xi,X2) represent detection non- uniformities across the detector acceptance m pt, ^ and 77. 
They average to zero, by definition, over the acceptance of the apparatus. For "reasonably behaved" detectors, one 
typically finds 

£2,0 - (£1,0)^ (56) 

Let us assume further that Si ^ 1. One then finds the raw two-particle momentum correlation is approximately equal 
to 



ptZt"' (Ar/, A^) « p^P^^P^- (Ary, A^) {1 + a2(A77, A^) - b2{ATj, A^)} (57) 



where 



. . J P2{Vl,V2,Vl,f2,Pl,P2)62{Vl,V2,fl,(P2,Pl,P2)'!^PlAp2dpidp2 

a2iAr],Aip) = — — — (58) 

J P2[Vl, 12,^1,^2, Pi, P2)dpidp2 

, , X J P2{Vl,V2,'Pl,'P2,Pl,P2)S2im,V2,'Pl,^2,Pl,P2)dpidp2 

&2(A7?,Av3 = — — — 59) 

J P2(Vl, 112,^1,^2, Pi, P2jdpidp2 

In general one has |ApiAp2| < 1, which implies 02 < 62, i.e. the difference 02 — 62 is non- vanishing. The inclusive 
transverse momentum correlation is indeed not perfectly robust as an observable. This implies a correction might be 
explicitly needed. A correction based on Eq. [55] is not practical because it requires knowledge of the two-particle 
density. It is thus simpler to carry an efficiency correction on a per particle pair basis, i.e. rather than counting "1" 
for each measured pair, increment the number of pairs according to 

^/£2{Vl,V2,(Pl,'P2,Pl,P2) (60) 

The above formulation of efficiency effects on the correlation neglects possible event detection biases incurred with 
small detection efficiency. Large particle losses may produce a significant modification of the pt spectrum of small 
multiplicity events. This is, thus, susceptible of significantly biasing the measurement of pt correlation oi p + p 
collisions. 

Figure [T] discussed in Sect. [V] displays PYTHIA events obtained with perfect efficiency and serves as a reference 
for simulations presented in Figures [7] through [5] obtained with limited detection efficiency. The effect of detection 
efRciency is obtained by individually and randomly rejecting particles on the basis of a given efficiency prescription, 

£2 = f{Vl,V2,fl,f2,Pl,P2) (61) 

One uses a two-particle efficiency equal to the product of single particle efficiencies, as per Eq. [55] The single 
particle efficiency is defined on the basis of a Fourier decomposition involving twelve terms as follows: 

e(<^) = eo 1 + ^ £i cos(n.^) (62) 



The coefficients are chosen to approximately model the azimuthal efficiency dependency of large detectors such 
as the STAR TPC and ALICE TPC. Spike structures are introduced in azimuth to simulate the effect of detector 
boundaries. Twelve sectors are used. The amplitude of the efficiency varies across the acceptance by ±10%. Figure 
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FIG. 6: Efficiency non-uniformities are sliown as a function of azimutfial angle, ifi. See text for details. 



ini shows efficiency non-uniformities introduced in azimuth. Figures [7] shows the robustness of the observables studied 
on the basis of Eq. [HH Our study of the impact of detection efficiency is based on single p + p collisions generated 
with PYTHIA. 

Figures [7]^ a-c) show correlations obtained with Eq. [51] using Eq ~ 0.8. For each produced particle one generates 
a random number, < r < 1. The particle is included in the analysis provided r < e. In Figures [TJ^a-c) we show 
the difference between the correlation functions obtained for perfect efficiency, i.e. e = 1 and reduced efficiency, i.e. 
(e) = 0.8, for all three observables of interest. The differences in the case of inclusive and event-wise observables are 
of the order of 0.5% only, whereas the difference observed for C is 1%. However, despite extremely small dependence 
on efficiency, we observe finite near-side and away-side structures in the difference. In contrast, C reflects the effect 
of detector boundaries. There are exactly twelve structures in azimuth which correspond to detector boundaries 
introduced as non- uniformities in azimuth. The differences are however numerically small and amount to ~1% of 
the signal at Aip = Arj = 0. We therefore conclude that all the three observables exhibit very little dependence on 
detection efficiency. 

We next include, for illustrative purposes, a linear efficiency dependency on the particle transverse momentum. 

e{ip,p) = £o (1 - ap) ^1 + X] ^' cos{nip)^ (63) 

Figures [8ja-c) show the difference between the correlation functions obtained for perfect efficiency, i.e. e = I and 
efficiency dependent on particle transverse momentum as per Eq. [S3] obtained with Eq = 0.8. The dominant features 
observed in Figures [S]Ja-c) remain the same as those observed in Figures [Tfa-c). We verified the robustness of these 
observables by further reducing the efficiency (Eq. [55]) to 50%, i.e. Eq = 0.5. The difference, pf^l^^^ {Ari,A(p) — 
P^e^^b (^Vj is of the order of 5% for inclusive and event-wise observables where as the difference for C increases 
by 2% only. We conclude that all three observables, inclusive, event-wise and C are quite robust. Their dependence 
on efficiency, even when it varies across the acceptance by as much as by 10%, is very small. 



IX. SUMMARY AND CONCLUSIONS 

We introduced formal definitions and notations for inclusive, event-wise and C correlation functions. We considered 
more specifically transverse momentum two-particle correlations expressed as functions of the particle pseudorapidity 
and azimuthal angle differences. We also studied a generalization of the observable C proposed by Gavin [l|, which 
he uses to determine the viscosity per unit of entropy of the nuclear matter produced in A 4- A collisions at RHIC 
and the LHC. 

We used the event generator PYTHIA [l^ to study and predict the shape of the correlation functions for minimum 
bias p-\-p collisions at y/s = 200 GeV. Simulations were carried out for particles in the ranges \r]\ < 1 and 0.2 < pt < 2.0 
GeV/c. The correlation functions exhibit a relatively narrow peak at A77 = Aip = 0, and a wide away-side band or 
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FIG. 7: (Color online) Difference of the inclusive, event-wise and C transverse momentum correlations function obtained with 
p + p collisions at -y/s = 200 GeV generated with PYTHIA for perfect efficiency, e = 1 and reduced efficiency, e = 0.8. (a) 
inclusive average p^^jL'^™ (A?7, Aip) — p^gL'o'? ('^'7> (b) event-wise average p^^L^^^ i^V; ~ pf^eLas^ ('^'7i ^v) a-nd (c) 
observable C, Ce=i — Ce=o.s- 



ridge extending over the full range of pseudorapidity \Ar]\ < 2. These features are qualitatively similar to those found 
in particle density correlations calculated with PYTHIA [25j . 

While the three observables are meant to determine the level of transverse momentum correlations between produced 
particles, their definitions differ significantly. One expects in the large multiplicity limit, values obtained with the 
three observables should be identical. However, in practice, multiplicities in a given bin of relative pseudorapidity 
and azimuthal angles {Ar],Aip) are rather small. This implies differences between these three observables can be 
substantial. We verified this expectation on the basis oi p -\- p collision simulations calculated with PYTHIA event 
generator [Tsj . We observed inclusive and event-wise observables to be approximately identical. Both, however, differ 
significantly from C. We conclude that for a precise comparative study of transverse momentum correlation functions 
for different colliding systems, at different energies and from different experiments, it is essential the same observable 
be used to characterize the correlation shape and strength. 

We studied the experimental robustness of the three observables. We studied the impact of detection efficiency 
using PYTHIA simulated p + p collision events with an arbitrary parameterization of the efficiency dependence on 
the detection angle and transverse momentum. We found that these observables are robust even when the detection 
efficiency reduces to 50%. 

We studied the scaling of the inclusive transverse momentum correlation function with the number of participating 
nucleons in A+A collisions assuming the p+p interactions are independent and with no rescattering of the secondaries. 
We found that, similar to other correlation observables the correlation strength exhibits a 1/Npart dependence for 
varying A -\- A collision centralities. 
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FIG. 8: (Color online) Difference of the inclusive, event-wise and C transverse momentum correlations function obtained with 
p + p collisions at -y/s — 200 GeV generated with PYTHIA for perfect efficiency, e = 1 and efficiency dependent on particle 
transverse momentum as per Eq. 1631 using eo = 0.8 and a = 0.05 (see Eq. I63p . (a) Inclusive average p^JL'^''^ (A?7, Aip) — 
p^gLas^ (At?, Avs), (b) event-wise average pf'^l^^^ {Arj,A(p) - pf'^l^^g {Arj,A(p) and (c) observable C, Ce=i - C'e=o.8- 

We additionally studied the effect of radial flow in p -\- p and A -\- A collisions on the magnitude and azimuthal 
dependence of these observables. We found that radial flow is responsible for the shift of the away side structure to 
the near side ridge like structure. 
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X. APPENDIX 



Experimentally, one measures at a given pt, rj, and ip a number of particles N{pt, rj, (p). This number is divided by 
the bin size to obtain the uncorrected density pexpiPt,Vi v)- general, the detection efficiency is a function of all 
three coordinates herein noted e{pt, fj, tp). The actual density, pTHiPt,VTf)j is obtained as follows: 

PTH {pt,V,'P) = — 7 ^ (64) 

The efficiency changes as detector conditions evolve over time. It may also be a complicated function of the measured 
event structure. It is the case for instance with a TPC in heavy ion collisions measurements. One observes experimen- 
tally that the detection efficiency depends on the detector occupancy, or particle track multiplicity in the detector. At 
low multiplicity, tracks are on average well separated and hence easy to identify and reconstruct. However, for large 
multiplicity, the distance between tracks is reduced such that they may cross and partially overlap. The efficiency is 
consequently a function of the event multiplicity as well as the beam luminosity. The dependence on luminosity may 
be corrected for by estimating e{pt, rj, ip) based on embedding techniques. We parameterize this dependence on global 
event factors with an index m. For convenience, one writes the detection efficiency at given multiplicity as follows: 

s(j)t,ri,ip\m) = eo{ri,ip\m)B{pt,ri,ip\m) (65) 

where eQ{pt,r],ip\m) represents the pt averaged efficiency, and B is a normalized response function. 



B{pf,ri,Lp\m)dpt = l (66) 
The number of particles detected at given rj and Lp is written as 



Nipt^V^f) = ^£o{v,'P\nT')B{pt,i],ip\m,)pTH{Pt,V,'P\m) 

m 



(67) 
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The measured average particle pt is then 

E ^oiv, 'P\m) J B{pt,T], (p\m)pTH{pt,V, v\m)ptdpt 
(Pt>exp (ri, = 1^ e^{r^, ^\m) J B{pt,v, 'p\m)pTHiPu V, 'P\m)dp.^ ^^^^ 



and is clearly not a robust quantity. However, in cases where the average pt is measured over a small range of 
multiplicity, one has 

/ \ f J B{pt,i],f)pTH{pt,V,'P)Ptdpt 

"'^P J B{pt, 77, ip)pTH{pt,V, fjdpt 

while this expression indicates the measured average pt still depends on the pt detector response, it is easy to verify 
that in cases where the pt dependence is mild, the above expression is a reasonable approximation of the actual 
average momentum. Consider for instance a case where the density has an exponential dependence on the momentum 
p oc exp(— pt/T) where T is a slope parameter describing the particle distribution. We evaluate the impact of the 
detector response using a linear approximation B (x 1 + apt- Assume T has a value of 0.6. The average pt changes by 
1% relative to the actual average pt when the value of the response parameter a is chosen to be 0.05. Smaller values 
of a lead to smaller deviations. One thus finds that unless a is exceedingly large, the detector response has a rather 
limited impact, a few percent only, on the measured average momentum. It is, nonetheless, a function of the detection 
coordinates rj and tp: deviations of measured (pt) may depend on 77 and </>. The above argument is straightforwardly 
repeated for measurements of correlations (ptPt) i^V^^v)- Given this quantity depends essentially on the square of 
the detector response, one expects the impact of a pt dependence should be rather small - although potentially a 
function of the particle coordinates iji, (pi, r/2 and ip2 ■ 



